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Abstract 

We evaluate the propagator of scalar and spinor in three dimensional quantum electrodynamics 
with the use of Ward-Identity for soft-photon emission vertex. We work well in position space to 
treat infrared divergences in our model. Exponentiation of one-photon matrix element yields a 
full propagator in position space. It has a simple form as free propagator multiplied by quantum 
correction. And it shows a new type of mass singularity.But this is not an integrable function so 
that analysis in momentum space is not easy.Term by term integral converges and they have a 
logarithmic singularity associated with renormalized mass in perturbation theory. Renormalization 
constant vanishes for weak coupling, which suggests confinement of charged particle. There exsists a 
critical coupling constant above which the vacuum expectation value of pair condensation is finite. 
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I. INTRODUCTION 



In three dimensional gauge theory infrared divergences is severer than that in four di- 
mension. In 1981,Jackiw first demonstrated it in the massless fermionic QED2+1 and after 
that other authors introduced the infrared counter terms to renormalize the infrared di- 
vergences [1,2]. There had been an attempts to solve Dyson-Schwinger equations to exam- 
ine the dynamical mass generation or chiral symmetry breakings in this model to improve 
low-energy behaviour [3,4]. Another important feature of this model is that it allows parity 
violating Chern-Simons term in the Lagrangean without violating gauge invariance[5].But it 
has not yet been clear the quantum effects of Chern-Simons term. During the same period 
infrared behaviour of the propagator in the presence of dynamical mass has been analysed to 
search the physical cut associated with massless photon [6, 7]. To maintain gauge covariance 
Delbourgo and Jackiw addded the vertex correction with Ward-Identities for soft particles 
in the determination of the infrared behaviour of the propagator with bare mass [9, 10]. In 
their works relation between bare and renormalized masses are given by the spectral func- 
tion. Their works are fundamental and important to understand the effects of soft-photons 
in comparison with perturbation theory. Of course their works includes non-perurbative ef- 
fects. In 1992 Delbourgo applied his method(gauge technique) to QED 2+ i and shown that 
massless loop correction to photon soften infrared divergences as TAppelquist et.al[3,6].In 



2 



reference [9] infrared behaviour of the scalar propagator in the presence of massless fields 
(photon, graviton) in four dimension was determined by solving the spectral function. In this 
work we apply the same method in three dimension to determine the infrared behaviour of 
the propagator. In section I we analyse scalar QED 3 and evaluate the scalar propagator in po- 
sition space based on spectral representation. In this case an approximaiton is made to choose 
the one-photon,one-meson intermediate state to derive the matrix element (Q\<f>(x)\n) .Most 
singular infrared part is assumed to be the contribution of soft photon emitted from external 
lines. Therefore we introduce photon mass to avoid infrared divergences. Lowest order spectral 
function contains mass and wave function renormalization up to logarithm in position space. 
Exponentiation of the lowest-order spectral function yields the full propagator in position 
space. It shows us a new type of mass singurality. Especially mass renormalization make a 
drastic change of the propagator. Its has a simple form as free propagator (with renormalized 
mass) multiplied by quantum correction. However this function is not integrable and it is dif- 
ficult to make a fourier transform of it. In the perturbative analysises the propagator have 
linear and logarithmic infrared divergences near on-shell.We mention the gauge transforma- 
tion property and see our solution satisfies Landau-Kharatonikov transformation. As far as 
the renormalization constant is concerned we evaluate it by position space propagator. The 
result is Z = for weak coupling,which is a signpost of confining phase in our approxima- 
tion. In section II we study QED 3 with spinor and see the spectral function for fermion is the 
same in section I. In this case there is a interesting possibility of pair condensation. Evaluating 
the vacuum expectation value (ipip) and we find that there exists a critical coupling constant 
(e 2 /(87rm) = 1) above which (ip-ip) remains finite. 

II. SCALAR QED 

First we assume the Kallen-Lehmann spectral representation of the propagator for mas- 
sive boson[l]: 
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A W = /^«p(*. I )< n |W(0)| n >, (1) 
a(s)ds 



Mir) . , 

J p l — s + le 

Im A{p) = a(p 2 ) = (2vr) 2 £ 5(p - p N ) (Sl\<j>\N) <iV|0+|fi) , (2) 

N 

|jV>=|r;A;i, k n >,r 2 = m 2 . (3) 

The spectral function is formally written as a sum of multi-photon intermediate states: 

-(p 2 ) = / ^E^(/^«(*») <w (**)Ew(i'--E*') 

x (fi|0|r i; fcx, (r; fci, ../c„|0 + |fi) . (4) 
Here we use the notation 

(/(*))o = 1, 

n 

(/(*))» =n ^ 

To evaluate the contribution of the soft-photons, we consider from the beginnings when only 
the nth photon is soft. Define the following matrix element 

T n (r; fci, .., k n ) = (ft|0|r; k 1: .., fc n ) . (6) 

We consider T n for k 2 ^ 0,we continue off the photon mass-shell by Lehmann-Symanzik- 
Zimmermann (LSZ) formula: 

T n = e^T n/i , (7) 
T% = J d 3 xexp(ik n -x)D x (Q\T<j)An(x)\r;k u ..,k n ) 

= J d 3 xexp(ik n ■ x) {n\T<f>f{x)\r; k u k n ^) , (8) 
where the electromagnetic current is 

fix) = le^^X^x) - 2e%0 + 0, (9) 
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here we assume the second term does not contribute in the infared behaviour of the propa- 
gator and the usual commutation relation for the first term 



8(x -yo)\jo(x),<f>(y)} = etf>{x). (10) 

Thus we omitt the second term in the present calculus. From the definition (8),T£ is seen to 
satisfy the Ward-Identity 

kniJ,T% — eT„_i(r; ki, .., k n -i). (11) 
Here we mention the simple proof: 

W = J d 3 xk n ^ex P (tk n -x) (Qim^lr-h,...,^) 
= i J d 3 xe^exp(ik n ■ x)d >1 (Q\T<j>f{x)\r\ h, 

= ee£{n\<j>\r;k 1 ,..,k n - 1 ), (12) 



provided 



d lt T{<f>{y)f{x)) = 5{x-y)e<f>{v), (13) 



where e is the charge carried by meson. To get the solution of the Ward-Identity we 
seek the non-singular contribution of photon for T n (r; k\, .., k n ). There is a possible solution 
which satisfy Ward-Identity [9] 

T^r-M-K) = ^+^ T w - 1 (r > fc 1 ...fc w _ 1 ),r 2 = m 2 . (14) 

Here we assume most singular contributions of photons emitted from external lines as in four 
dimension. These arise from Feynman diagram in which the soft-photon line with momentum 
kfi and the incoming meson line can be separated from the remainder of the diagram by 
cutting a single meson line. Detailed dicussion and evaluation of the explicit form of T n are 
given in ref [9]. Of course for n — 1 

_ e(2r + k) ■ 6 
ll ~ (2r + {b) 



and we can determine T n recursively with (14). In this case T n becomes 

e(2r + k) ■ > 
2r ■ k + k 2 



T. - I ) • (16) 



n 
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The one photon matrix element by LSZ is 

Ti = (in\T((f) in (x),ie / Sy^ n dv<f) in {y)Af n (y))\r; k in 



(2r + k),eHk) 



exp(i(r + k) -x). (17) 



(r + k) 2 — m 2 

From this lowest order matrix element the function in (4) is given 



F= (^l0(^)k;A;)<r;A;|0 + (O)|fi> 

one photon 

and a is expressed by 

d 3 



/d x f d t 

j^-^ exp(ip ■ x) J — exp(ir -x)exp(F). (19) 

Here is a polarization sum we have 

n„ = -(^-^)-d^, (20) 

and the free photon propagator 

^ = ^#*"-^ + ^l- < 21 > 

we get 

d 3 A; 



/d 6 k 
— eM^-x)6(k ) 

(^r+krm s(k 2 ) 

X [ (2r-k + k 2 ) 2 +[d L) k 2 1 



(22) 



It is natural to set 5(k 2 )/k 2 equals to — 5'(k 2 ).We see the reason in calculating the one loop 
self energy of meson diagram, and take imaginary part. This expression is the same as the one 
given above with the interpretation S(k 2 )/k 2 = —5'(k 2 ).We must introduce a small photon 
mass /i as an infrared cut off. Therefore ( 22 ) becomes 

d 3 k „.,„. ,., WI - /r9 , w m 2 1 , . , _ . d 



-e 2 



I (2^° )eXp( ^^^^^ (23) 



The evaluation of the integral to F is described in Appendix. The result is 

F= ^ d ' 2) + ^- + &r to <"*> - ^ - i x(i + -<~ 2) ' (24) 



It is helpful to use position space that it shows us easily the short and long-distance behav- 
ior. The relation between perturbative spectral function and F(x) is 

exp{ip ■ x) J — exp{ir ■ x)F(x). (25) 

Here we think about the physical meanings of each term. 

exp(— mx) 

cr{x) = — exp(F(x, e, //, m)), (26) 

a(p) = F.T.(a(x)) = J -0^ exp(ip • x)a(x). (27) 

The terms proportional to x correspond to mass renormalization and others to wave 
function renormalization. 

^)=^(A f M - { Z x +B)X) (^- CX+D ' < 28 > 

(d-2) + ^—,B = ^-(d + 1 -2),C = |-, D = -— . (29) 



87TyU 87rm 87r 87r 87rm 

From the form of cr(a;),wave function renormalization effects denoted by constant and 
D\n.{fix) modifies the dimension of the propagator as anomalous dimension in four dimen- 
sion. Mass renormalization is proportional to x and x ln(/xx). Logarithmic corrections descrive 
the logarithmic infrared divergences and adds a drastic effect as a factor {jjt,x)~ Cx .This factor 
(/ix)~ Cx is a new type of mass singularity in three dimension. In perturbative analysis these 
will be shown clear. By the definition (1),(2),(3) a(x) is a full puropagator.Free propagator 
in three dimension is 

, ^ x exp(— (m + B)x) . 1 . . 

a {x, m + B) = U-, ao p, m + B )= 30 

4nx (m + B) z + p z 

and its quantum correction is expressed by exp(A)(fj,x)~ Cx+D .We find in this expression that 
a(x) is finite provided 

\a(x) = exp(A) eXp( ~ ( 4 ^ + B)X) ^x)- Cx+D \ = finite, (31) 

POO 

< / a(x)dx < M, < D, (32) 
and there exists a(p) 

a{p) = r x -^i a(x)dx . (33) 

Jo V 



7 



In this case we have not a simple pole even in the Yennie gauge d = 2 in which we see the 
singularity at p 2 = (m + e 2 7 /87r) 2 .Next we treat the momentum dependece of the propagator 
a(p) (32). The integral is not analytic and we cannot get the precise expression for a(p). 
If we expand (fj,x)~ Cx+D 

{fix)- Cx+D = 1 + {-Cx + D) ln(//x) + \{-Cx + D) 2 (\n(iix)) 2 + \{~Cx + Df{\a{^x)f + .., 

(34) 

term by term integral of a(p) converges. Using the following integrals 

h — f — - exp(— mx) ln( /ix)dx 
Jo V 



o V 

7 ln((m 2 + p 2 )/fJ. 2 ) ln((m — \/—p 2 )/{m + a/— p 2 )) 



p2 _|_ m 2 2{p 2 + m 2 ) p 2 + m 2 



(35) 



I 2 — [ exp(— mx)x In (fix)dx 

Jo P 

[l n((m _ y/q?)/(rn + v^)) + \n((p 2 + m 2 )//i 2 ) - 2(1 - 7 )], (36) 



(p 2 + m 2 ) 2 
a{p) up to 0(e 2 ) is given 

A mB 



°" (2) (P) = [-T— 2 + 



m 2 + p 2 p 2 + m 2 2(p 2 + m 2 y 



+ Dh-CI 2 \,p=yf^P. (37) 
In this case the renormalization constant Z (residue of the pole term) becomes 



Z = 1 + A + L>(- 7 + - ln(^— )) p2 _ oc - oo. (38) 



In the Yennie gauge there remains a renormalization 



e 2 

Z = l + — ln((p 2 + m 2 )/// 2 )^^ -> oo. (39) 

lo7rm 



In Minkowski space, p 2 — ■> — p 2 ,with dicontinuity 



— !— = P.y- + i7r5(x) (40) 
x — le x 

( 1 ., 2 =P.V- 2 +™5\x), (41) 
(x — «e)^ x" 1 

— — ln(x - ie) = inS(x) ln(x) + P.V-m. (42) 

37 37 
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we can determine the the strucure near p 2 = —m 2 by the imaginary part of a(p).Here we 
notice the second-order spectral function (p) 



^ = exp( mx) +A _ Bx + {D _ Cx) ^ ))? 



(43) 



(2) 1 , A- 1 D 2mB 



m? + p 2 m? + p 2 (m 2 + p 2 ) 2 
D H^{m 2 +p 2 )/^ 2 ) 
{p 2 + m 2 ) 

+ , 2 ^ 2 , 2 N(P 2 + ™W) - 2(1 - 7 )]. (44) 

Here we look at mass renormalization part 

_ m(2B + 2C(1 - 7) - C ln((p 2 + m 2 )//i 2 )) 

m 9 i 9 

^2 _|_ m 2 

_ e 2 2(rf-l)-ln((p 2 + m 2 )/^ 2 ) 
87r p 2 + m 2 

Spectral function lma^(—p 2 ) reads 



(45) 



T (2) / 2 \ 2 

[ ~ P ' = {l + A- 7j D)5( P 2 - m 2 ) - ^-m{d - 1)5' (p 2 - m 2 ) 

7T 47T 

D6(p 2 - m 2 ) Cm6(p 2 - m 2 ) 
m 2 — p 2 (m 2 — p 2 ) 2 

Usually we estimate the second-order spectral function in the Feynman gauge d — 1; 
W2)(V) = [(l-^)^ 2 -m 2 ) 

7T d=l ^ 8lTfJ, J KF ' 



^2 _ m 2 ^ m 2 _ p2\2 



This shows the renormalization constant Z has an linear infrared divergence as \i — > O.If 
we add higher order corrections linear infrared divergences cancell each other in the same 
mechanism in four dimension[15,16,17]. In the Yennie gauge it has cuts 



Ima^i-p 2 ) 2 2 e2m x'f 2 2n 
— oyp — m ) — o (p — m ) 

IX d=2 47T 

+ V-(^—2 + { 2 m \^ V 2 ~ ™ 2 )- ( 48 ) 
8nm p z — m z (p 2 — m 1 ) 1 

To see the singularity structure near p 2 = — m* 2 it is better to expand around p 2 = — m* 2 .We 
get 

a {2 \p) = [ — 1 — - + (expM) - 1) x {-CUm -> m*) + Dh{m -> m*)]. (49) 

m *2 _|_ ^2 



In this case the renormalization constant Z : 

Z = l + (exp(A) - 1) + D(- 7 + i ln(^±^)) (50) 

is divergent at p 2 = —m* 2 . There is an interesting contribution from exp(F 2 ) 
sin (pa;) 



Jo P 

cos(7ra 



exp(— mx) (fxx) a dx 



= -T(a + l)i^ipV + m 2 )- 1 -"' 2 ^ 



\/—p 2 \/ —P 2 + m \[—p 2 — m 

2 



~ i\/ —P 2 — m) 1 a near p 2 = — m 2 , a = . (51) 

Here we notice that this type of singularity also appears in four dimension where a = 
—a(d — 3)/(27r) but in this case a is gauge invariant. Formally fourier transform of a{p) can 
be written by double fourier transform 

<P) = J ^ F - T - C M ~ mX) W )(g) x F.T.(»x)- c *(p - q), (52) 

but we do not discuss details of it in this paper. As we see in the annalysis in terms of e 2 ,the 
contribution of the function (fix)~ Cx leads to a new type of mass singularity ;divergent series 
at p 2 = — m* 2 in o"( 2 )(p).We have seen that the fourier transformation of the propagator is 
very difficult.lt is interesting to study the phase structure of the model. First we estimate 
renormalization constant by the following sum rule[10,ll], 

Z- 1 = / a(u)du>, (53) 



ry-i ,. 2 f , n, sm(px) exp(-ujx) , . A . 

Z = hm p 2 / a(uj)dw — , L dx, (54) 

p 2 ^oo J p Anx v ' 

lim F(p) = lim r nXSm ^ F(x)dx = lim ^-(xF(x))(0), (55) 

p^oo p— >oo Jq p p— >oo 

we can evaluate this quantity by direct substitution of a(x) into the above equations and 
take the limit 

Z^ 1 = exp(A) lim exp(— mx){^x)^ Cx+D 



0. (56) 
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Next we consider the gauge dependence of the propagator[12,13,14].When we write the 
photon propagator as 

D^(k) = D$(k) + k^KM(k), (57) 
we seek the change of the propagator. Under the gauge transformation defined: 

(f>(x) — > exp(iex(x)), (fi + (x) — > (p + (x) exp(—iex(x)), 
A^x) - A^x) - d, X , 5{<P + X<P) = -d,x^ + X<P ~ 2^0+0^), (58) 

propagator changes 

2V - i {T[A,{x) - dZx(x)][Mv) - dliv)]) 

= D^ + i(Td?x(x)dy(y)), (59) 
5D, V = dffiMix -y),M(x-y)=i (T X (x) X (y)) 
A -> -% (T(f)(x) exp(iex(x) exp(-iex(y)(f> + (y)) 
=A (T exp(iex(x)) exp(-iex(y))) . (60) 

Here we expand in 5x 

(T exp(ie5x(x)) exp(-ieSx(y))) = ~^e 2 (5x 2 (x) + 5x 2 (y) - 2T5x(x)6x{y)) , 

5 A = ie 2 A(°\x - y)[M(0) - M(x - y)\, 
i[M(0) - M(x -y)] = - (T5x(x)5 X (y)) • (61) 

Usually M is identified as the gauge fixing term 

M(k) = -d/k\ 
. 9 , , . . . 9 , f d 3 k exr>(ik • x) — 1 

= " aF* < 62 > 

with subtraction of infrared divergence at k — 0. We get the propagator(59) in covariant 
gauge for three dimension 

e 2 d 

A(d, x) = exp(-— x) A(0, x) (63) 
e 2 (i 

= exp(- — x)a(0,x), (64) 
11 



Since cr(x) is a full propagator it obeys the Landau-Kharatonikov transformation in our 
approximation except for linear infrared divergent factor A 

. , exp(-(m* + de 2 /87r)x) \-Cx+d ,, r > 

a(d,x)= —exp(A)(jix) ° . (65) 

Therefore we notice that the Landau-Karatonikov transformation is regularization depen- 
dent. We keep the Landau gauge and not take the covariant d gauge to avoid longitudinal 
photon. Here we mention the radiative correction of mass to understand the mass shift. In the 
Landau gauge our approximation changes the mass mtom* = m + e 2 (7 — 2)/(87r).The con- 
stant 7 appeared after the expansion of Ei(l,/ix) in //.In ref[5],gauge dependent self-energy 
is estimated at 0(e 2 ) 

E(m)=d|^, (66) 
and the physical mass for fermion defined in the Landau gauge 

m phy = m + S(m). (67) 

If we estimate one-loop self-energy in the Landau gauge 

S(m) = i-, (68) 

which is two-times larger than e 2 /(4ir) in m* since we used D + (x) instead of D(x).We see 
the same effects occured for scalar with vertex correction. We call m* = m + e 2 (7 — 2)/87r 
the renormalized mass in our approximation. At first sight, mass shift seems to be peculiar 
but infrared behaviour is governed by one-particle intermediate state. 



III. SPINOR QED 

In this section we study the spectral function for massive fermion. Similar to scalar case 
propagator and the spectral function are defined 



Sip) = I d 3 xexp(ip ■ x) (n\T(^(x)ilj(0))\n) 
\ p{uj)duj 



(69) 

p ■ 7 — uj + te 

p(p 2 ) = (2vr) 2 J2S (3) ( P -Pn) (N\m). (70) 



N 
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Matrix element is 



T n = (n\^\r;k u ..,k n ), (71) 
T n=~ J ^x^{ik n -x){Q\T-il;f\r-M,...K-i), (72) 

provided 

U x T^{x) = TijD x A,(x) = T^(-j^x) + d%d ■ A{x))). (73) 
In the similar way to the scalar case T n satisfies Ward-Identity: 

dlT^j^x)) = -e^x), (74) 
^T(^(x)) = e^(x), (75) 
k n ^(r, k u k 2 , ..k n ) = eT n _i(r, k ± , k 2 , ..fc n -i), r 2 = m 2 . (76) 

One photon matrix element by LSZ is 

Ti = ( in\T(ijj in (x),ie / d 3 yip in (y)-f fl ilj in (y)A^ n (y))\r; k in 



le 



J d 3 yd 3 zS F (x - y)7 M 5 (3) (y - z) exp(i(k -y + r- z))^(k, X)U(r, s) 



= -^ «\ k L 1 + Z ^(k, A) exp(.((A; + r) • x))U(r, s), (77) 
{{r + k) A — m z ) 

where U (r, s) is a two-component free particle spinor with positive energy: 

^U(r, 8 )U(r,8) = l^^. (78) 

In this case the function F becomes 

2 f d 3 k \nn0un2u_r( r "'"^)'T + m 7 • r + m (r + A;) • 7 + m 



(79) 



F = ?t(-- + < l - 1) - x ln <^» + + -r) + f(d-i)(--x). (80) 

»7rr z /x 87rr 87r /1 

Thus the lowest order spectral function is exactly the same with that in the scalar case 
except for the normalization factor in the phase space integral: 
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exp(ip ■ x)p(x) 

/d^x f in 

- — — exp(ip-x) / d 2 r — exp(ir • x) exp(F). (81) 
[2ix) 6 J r 

After angular integral we get 

p(p) = -^exp(A) f°° dx sin(px) ex ^ B ± m ^ (^)~ (82) 

Therefore the structure is the same with scalar case as we mentioned in the last section. 
Consequently there is no infrared divergences at p 2 = m 2 . 

For the renormalization constants we apply the same argument as for the scalar,full 
propagator is expressed in the dispersion integral 

/GXJ) f — UJOC J 
p(uj)duj(i^ ■ d + uj) = — (ry • daix) + cr(x)), 

Here we use the fourier transform of Sf to determine Z^.For free case the propagator in 
position space is 

S - (x) = (!7 ■ a) W^P = M^W = -7^5?' (83) 

We see the dimension of the Sp\x) is equal to l/x 2 .I momentum space we obtain 

r xsin ^ i dx = (84) 

Jo V x 2 2p' 
and this shows the ordinary expression in momentum space 

4 „, = 7_P = 7J.1 ( 



Z? = lim [(. 7 ■ da(V^))^=] = hm ( ^^-f ) (86) 



Trexp(A) lim j, ( ex P( ™r) Cr+D 
x^o+ dr r 



(1 < D) 
00 (1 > D) 



(87) 



rf.exp(-mr). x _ rr ,n, , w . r ,. +D[ m + C D — l C\n(pr) 1 
-r(—^ - H Cr+D = exp -mr /ir C?+D + ^ 

,r = v /Z a?. (88) 

moZz 1 = / ujp{u)duj = exp(A) lim (/ia;)~ Ca;+D = 0. (89) 
J x ~*°+ 
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Therefore there is a confining phase; Z 2 = 0(Z 2 _1 = oo) of charged particle for weak coupling 
constant. In this case bare mass vanishes for all coupling. Order parameter for the vacuum 
expectation value of pair condensate is given 



M) = -itrS F (x) = -2m 2 exp(A) lim (exp(-77ix)(x) -c,x+£> - 1 )) 



= -2m 2 



(1 < D) 
finite (1 = D) 
oo (1>D) 



(90) 



Here we notice that there is a critical coupling constant D cr = l(e 2 /(87rm) = 1). 



IV. SUMMARY 



We have seen how the Ward-Identity for soft photon may be applied to three-dimensional 
electro-dynamics to extract full propagator in position space and the infrared behaviour of 
the propagator in momentum space. Neglecting unconventional terms in the elecrtomagnetic 
current for scalar,spectral functions for scalar and spinor coincide each other. Exponentiation 
of the lowest order spectral function corresponds to the infinite ladder approximation. Since 
the theory is super renormalizable the mass are corrected to add some finite radiative cor- 
rection. We found a new type of mass singularity in three dimensional QED in position space 
as (fix)~ Cx . In momentum space it shows us logarithmic mass renormalization and yields a 
singular infrared structure of the propagator. Wave functuin renormalization D plays the 
role as anoumalous dimension .These are the consequences of logarithmic infrared diver- 
gences. Renormalization constant vanishes for spinor case and there is a confining phase for 
weak coupling. This picture is consistent with perturbative logarithmic infrared divergences 
at p 2 = m 2 .lf the coupling constant is smaller than D cr (D = e 2 /8irm < 1), vacuum expec- 
tation value (ipip) becomes finite. In our lowest order spectral function there remains linear 
infrared divergences that was regularized by photon mass.O(e 2 ) correction to the external 
line and the cancellation of infrared divergences are now in progress[17]. 
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VI. APPENDIX 

To evaluate F we use the function D + (x) in three dimension: 



DJx) = — ^- [ exp(ik ■ x)9(k°)5(k 2 - fi 2 )d 3 k 
(2n) 2 i J 



irkdk exp(-ax) 
Mkx)—== = P Q L: \ (91) 



{2nytJ uv '2y/k 2 + fj? Smx ' 
with the following parameter trick as in ref[9] 



lim / daexp(i(k + ie) ■ (x + or)) = exp ^ fe ' x \ (92) 
e->o+ y k-r 

lim / ada exr>(i(k + ie) ■ (x + ar)) = — ^— -r — . (93) 
6-o+ J KV (k-r) 2 

The function F is written in terms of the parameter integrals 

/»oo /»oo 

F = ie 2 m 2 / adaD + (x + ar, fi) — e 2 / daD + (x + ar, //) 
Jo ./o 

-ie 2 (d- l)—D + (x,n) 

= ^(- 6XP( ~^ ) +^Ei(l,^)) - ^-Ei(l,/xx) + (d- l)-^exp(-^), (94) 
87rr z /i 87rr 87171 

where the function Ei(n, x) is defined 

, N f°° exp(-xt) , , . 

Ei(n,x) = j PV tn > dt. (95) 

It is understood that all terms which vanish with ji — > are ignored. 

Ei(l, /ix) = -7 - ln(^) + O(^), (96) 
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p 2 m 2 1 

Fi = + x(l - ln(^) - 7)) + 0(/i), 

F 2 = ^-(ln(/ia;)+7) + 0(/i), 
F s = f (--x)(rf-l) + 0(/i). 



Here 7 is Euler's constant. Using the integrals 

J d 3 xexp(ip ■ x) J d 3 r5(r 2 — m 2 ) exp(ir ■ x)f(r) = f(m), 



/d 3 x exp(ip ■ x) f - — —5(r 2 — m 2 ) 
J (2vr) d 



m 2 + p 2 



we get F in position space 



F = -^—(d - 2) + ^— + ln(ax) - la(ax) - ^—x(d + 7-2). 
8tt/z ; 87rr 8vrr VP ; 8tt vp ; 8tt v ' y 



After integration over r 



exp(ir • xcos{6))d6 = nJ (rx), 

nrJ (rx) exp(-mx) 



I 

Jo 



dr- 



y/r 2 + m 2 x 
we obtain the full propagator 

, . , exp(-(m + B)x) N _rv+D 

<t(x) = exp(A) ^ K v —(fix) Cx+D } 

Attx 

^=^(d-2) + ^:,S = ^(d + 7 -2),C=^, J D ( 



87171 87rm 87T 87T 87rm 
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